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1NTRXXCTI0N 


Hie ability to measure ocean currents and certain ocean wave and wind speed character- 
istics over larqe ocean areas and long periods of time could have significant impict upon 
open ocean and coastal activities. Ship routing, search and rescue operations, oil .spill 
modelling, meteorological and exx\mographic research, and recreational activities are just 
a few of the areas whore quick and reliable information is desired. With the advent of 
satellite altimetry, it is now possible to estimate nvmy of these ocean characteristics 
with surprising accuracy. 

Three products which can bo estimated from satellite altimetry have particular impor- 
tance. They art' geostrophic ocean current petition and velocity, significant wave height 
(SVJH) # a statistical parameter which relates the height of the ocean raves to mean sea 
level, and surface wind speed. All of these products can bo conputod on a global basis and 
in near real-time. 

The National Aeronautics and Space Administration (NASA) Geodetic Experimental Ocean 
Satellite (GFOS-3) spacecraft ( fit jure 1) was launched from the Air Force Western Test 
Range on April 9, 197.8. Hie satellite orbits the earth every 101.8 minutes in a near- 
circular, 65-degree retrograde orbit. Hie primary instrument on board the 340 Kg satellite 
is the radar altimeter, developed for NASA by the General Electric Corporation. The altim- 
eter operates at a frequency of 13.9 GHz, transmitting 100 pulses per second. The pulses 
transmitted by the altimeter reflect from t he earth's surface and are received by the 
spacecraft yielding an altitude measurement with an RS IS accuracy of 70 cm. This configura- 
tion gives a footprint size of 3.6 km wide and 11 km along the track from an orbit of 
840 km with zero eccentricity. A detailed technical discussion of the GEOS- 3 radar altim- 
eter is given by tofmoistor et al . , (ref. 1). 

The altimeter is instrument al with 16 sample and told gates which prxvide information 
about the shape and amplitude of the return waveform. Hiis information can be used to 
determine a number of interesting and use fill parameters including water- land and water-ice 
boundaries, ocean surface wind speed , and significant: wave height. At the insistence of 
the 01X6-3 Project Scientist, H. E. Stanley, an experimental version of a SWH estimating 
algorithm was built into the alt imeter processing software (ref. 2). Several refinements 
to the original SWH algorithm wore produced as a result of the analysis of early GEOS-3 
on-orbit data. Hie analysis also led to the development of algorithms to provide ocean 
sea surface wind estimates and the water-ice a to water- land boundary estimates* 

These estimates proved to to so useful that, in 1978, NASA established the GEOS-3 
Near- Real-Time Data System for disseminat ,ug Gulf Stream boundary positions and world- 
wide significant rave height and wind stood estimates. In the near-real-time system, the 


GEOS- 3 data were acquired from the 13 NASA Space Tracking and Data Network locations or 
the ATS- 6 satellite and transmitted to the NASA Goddard Space Plight Center (GSFC) in 
Greenbelt, Maryland , where they were buffered to magnetic tape in real-time. The buffered 
data were then transmitted to the Computer Sciences Corporation INFONET center in 
Beltsville, Maryland, where significant wave height and wind speed were computed and made 
available to user-supplied terminals on a call-up basis. At the same time, the buffered 
magnetic tapes were sent to NASA V&llops Plight Center (WFC) , where, using predicted or- 
bits, an estimation of the position of the Gulf Stream was calculated. This information 
was then also sent to the INFONET center. 

Due to the efforts and the early recognition of the potential of the GEOS- 3 radar 
altimeter data by the Bureau of Land Management (BLM) Outer Continental Shelf Office and 
in particular by Dr. E. D. Wbod, an interagency agreement was reached between BLM and 
NASA. This report represents the first application of GEOS- 3 data for a large, comprehen- 
sive environmental study such as BLM's South Atlantic Bight Study. 

DATA PROCESSING FOR SEA SURFACE TOPOGRAPHY 


The study of the general ocean circulation has long occupied the central stage of 
physical oceanography and the most important phenomena in the ocean circulation investi- 
gations are the major boundary current systems such as the Gulf Stream and Kuroshio. Bor 
the Gulf Stream, most of the results from past studies have been summarized by Stommel 
(ref. 3) in a classic work. Although the results are voluminous, our knowledge of the 
global ocean circulation in general and the Gulf Stream in particular is still quite 
limited. Recently, the development of an accurate radar altimeter operational under all 
weather conditions has reached a new level of success. Direct measurement of the sea sur- 
face topography with an accuracy of 20 cm for one second averaging from a stable moving 
reference platform provided by a satellite is new possible. This offers a direct way to 
measure the dynamic topography which until new could only be inferred from hydrographic 
data. The measurements made by a satellite altimeter are illustrated by figure 2. An 
ideal altimeter will yield measurements on satellite altitude which may be expressed as: 
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where h is the satellite altitude measured by the altimeter, h is the satellite height 
a s 

above a reference spheroid determined from the satellite tracking data, h is the geoidal 

y 

deviation from the spheroid due to gravitational anomalies and Ah is the height deviation 
from the geo id. dr * to the dynamic processes in the ocean such as the currents and the 
tides. From the altimeter data, the surface current can be calculated directly from the 
surface slope using the geostrophic equation (ref. 4). 
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The overall data processing is summarized in figure 3 in the form of a flowchart. An 
attempt to more fully explain each step is presented here. The initial task consists of 
obtaining a search of the data files for passes within the study area. Many of the passes 
are eliminated due to a variety of reasons, among them: lack of orbital data, poor data, 

data obtained over land only, and data of only a few seconds duration. Once a set of data 
is initially obtained, each GEOS- 3 pass within that set is identified by orbit number, 
time, latitude and longitude. This information goes into the Dynamic Ocean Surface (DOS) 
computer program which is designed for the purpose of obtaining a sea surface profile from 
the altimeter height data. 


Sea Surface Height Profiles 

Averaged altimeter measurements are preprocessed and converted to sea surface heights 
as defined by Leitao et al. , (ref. 5) and treated here as raw sea surface heights. The 
preliminary data processing procedure used for this study by the DOS computer program is 
summarized in figure 4. The raw sea surface heights shown in A are first edited to elimi- 
nate anomalies due to internal instrument noise. The edit criterion is based on a predict- 
ed sea surface height calculated from fitting a straight line through the last eight sec- 
onds of data (80 data points representing 60 km in physical distance). Any point differing 
from the predicted height by more than two meters (approximately three standard deviations 
of the noise level of 70 cm) is replaced by the predicted height value. Results point out 
that less than 1% of the calculated predicted data are used. The edited sea surface height 
data shown in B are then filtered using an 81-point, equal weight, mid-point filter. This 
filter is chosen since it reduces the noise without seriously (compromising the sea surface 
signature. The eight seconds are selected so that the accuracy could be maintained below 
the 10 cm level; and the noise level is sufficiently reduced prior to differentiation. 

After the filtering process, the smoothed sea surface height is referenced to the Mader 
(ref. 6) 5 1 x 5* geoid (See C). Subtracting the geoid from the smoothed sea surface height 
results in a residual which is nearly flat in the open ocean as shown in D. Next, to mini- 
mize the error between the geoid and the smooth sea surface height south or east of the 
mean position of the Gulf Stream, a linear fit is made to the residuals over the section 
representing the open ocean. The straight line is then subtracted from all residuals thus 
removing any potential orbital bias or slope errors and producing an estimate of the sea 
surface profile as shown in E. Cnee all the sea surface profiles are obtained for each 
calendar month, they are compiled onto monthly tapes. Thus, thirty-six monthly tapes are 
available. 

The determination of where to fit the data to the geoid is first done blindly. A 
second fit is then perform'd once the first attempt is analyzed. Since temporal changes 
in the sea surface could have affected this second or a-posteriori fit, it is decided to 
further analyze this problem. At this point all passes within the study area are grouped 
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geographically by bands running parallel to the GEOS- 3 subsatellite tracks, each band 
having a width of 2,5° of latitude at the equator (figure 5). Each band is subsequently 
analyzed independently of the others and each individual pass within each band inspected 
to ascertain the maximum southern excursion of the Gulf Stream within each band. Once 
this analysis is performed, it is found that by dividing each band into five sub-bands the 
extreme excursions are better modelled. The difference in time and position between a 
determined parallel and the location of the maximum southern excursion of the Gulf Stream 
is calculated for each pass within a sub-band. The time for the termination of the fit to 
the geoid iy then calculated geometrically ( f Vjure 6). 

About 25% of the GEOS- 3 passes are octree ' -xi by tMs method and the new, corrected, 
monthly compiled data tapes are evaluated. s Fo facilitate the data handling, all data are 
stored on disk .it- a rate of every tenth measurement. A short study was conducted to ana- 
lyze the differences between the point to point sea surface profiles and the every tenth 
profile (figure 7). The results point out that there is no significant difference 
between the tenth second and the one second data. Hence, one second data are used. 


Sea Surface Topographic Maps 

With all the data stored on disk, each month's data are processed by first using the 
MODSEAHT computer program, a modification of SEAHT (ref. 7). MODSEAHT carbines all the 
monthly profiles using a minimum variance technique to minimize the differences in height 
at the intersecting points. The adjusted values of the topographic height which minimize 
the errors at the intersecting points are treated as the directly measured dynamic heights 
for the month. These adjusted data are then reticulated for automatic contouring. This 
technique divides the scanned area into 1/2 degree cells. All adjusted measurements 
occurring within a cell are averaged and cells with no observations are filled by averag- 
ing two-way linear interpolations, one along a constant band of latitude and the other 
along a constant band of longitude. Once all the cells are filled, blank cells are filled 
with a coded -9999, and an automatic contouring technique is applied. This produces sea 
surface topographic maps with contour lines from +200 cm to -200 cm with a contour interval 
of 20 cm. 

Using the 1/2° x 1/2° contour maps, a 1° x 1° averaging of the sea surface flow is 
performed. The geostrophic component of velocity is calculated and positioned according 
to the contours as observed in the topographic maps. The calculations are based on the 
simplified geostrophic equation: 
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2 sin c 
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where V is the geostrophic surface velocity; g is the gravity acceleration of 980 cm/sec ; 
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p, the density of sea water taken as 1 g/cm ; 8 the angular speed of the earth; <}>, the 
latitude; Ah, the total height difference; and AL, the horizontal distance lie tween contours. 
Graduated arrows are then placed on the map indicating the magnitude and the direction of 
the sea surface flow. 


DATA PROCESSING FOR SIGNIFICANT WAVE HEIGHT AND WIND SPEED 

Various experiments performed at WFC (refs. 8, 9, 10) and at other sites by other 
principal investigators (refs. 11, 12, 13) indicate that the GEOS- 3 data canpare favorably 
with other estimates and measurements of sea state and surface winds. The first large 
scale comparison of GEOS- 3 data with independently derived sea state information present- 
ed several variations of the c * iginal significant wave height algorithm and estimated the 
accuracy of the altimeter-derived significant wave height to be 55 cm (ref. 14). Present 
studies at WEC estimate the accuracy of the altimeter-derived wind speed to be 1 m/sec. 

The algorithms derived in this section have been in use at WFC in preprocessing the GEOS- 3 
altimeter data since duly of 1975. The same algorithms were used at GSFC in the near-real- 
time data network established in March of 1978. 

Significant wave Height Estimation Algorithm 

The model used to obtain the waveform geometry for negligible significant wave 
heights is based on the geometry of a square pulse impinging upon a flat sea surface as 
illustrated in figure 8. A pulse of duration T leaves the spacecraft (observer) at time 
t = 0, and traverses a distance H to the sea surface, arriving at time t = t^, where it is 
reflected back to the satellite at time t = t£ = 2t ^ . Therefore 


. 2H 
fc 2 = — 


(3) 


where c is the speed of light* At t = t^ = t 2 + t where t < T, the square pulse is 
observed impinging upon the sea surface, where an expanding circular area is illuminated. 
The radius r of the area is related to H by 


r 2 = (H+h) 2 - H 2 = 2Hh + h 2 « 2Hh 


(4) 


since n << 2Hh. The time t = t^ is given by 


t = t + T = ^J+T 
•J 0. C 


(5) 
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aid corresponds to the two-way travel time between the satellite and point P on the 
surface. Therefore, 


35 


,2H , . 

(— + T ) 


H + h 
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or 



h 
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which yields 



Combining equations ( 4 ) and (8) yields 


r 2 = Hct 

and the area of the illuminated region is given by 

2 

a - irr = itHct 


( 6 ) 


( 7 ) 


( 8 ) 


( 9 ) 


( 10 ) 


Note that the observed area of the illuminated region, or equivalently the power received 
by the satellite, increases linearly with time. When t becomes greater than t^, the 
trailing edge of the pulse has reached the sea surface and the area illuminated becomes 
an annulus with inner radius and outer radius g given by 

r 5 = Hctj. (11) 


r 


2 

6 


= Hct- 
6 


( 12 ) 


and width W and area A given by 


W = r 6 - r 5 = Hdt 6 - Hct^ (13) 

A = nr 2 - irr 2 = irHc(tg - tg) = irBcT (14) 

Since the area of the annulus remains constant, the power received at the satellite re- 
mains constant until the antenna-beamwidth effects cause the power received to decay. 
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The characteristics of a square pulse impinging upon a sea surface with negligible 
SWH as seen from the satellite can be suimrarized in four stages; (1) no power is received 
until the leading edge of the pulse is observed striking the sea surface, (2) after the 
leading edge of the pulse reaches the surface and before the trailing edge does, the power 
received increases linearly with time, (3) after the trailing edge of the pulse reaches 
the surface, the power received remains at a constant plateau value, and(4) after the 
antenna-beamwidth effects become non-negligible, the power received begins to decay. 

These four stages are depicted in figure 9, which represents the idealized mean return 
pulse shape. Due to the nature of the altimeter, the instantaneous power received fluc- 
tuates, making it necessary to average a number of pulses in order to determine the mean 
pulse shape. The GEOS-3 satellite receives the return pulses in 16 waveform sampling 
gates. These 16 values, called Instantaneous Return Samples (IRS), are collected 100 
times per second by the satellite and averaged on-board in an attempt to construct an 
accurate representation of the mean pulse shape. The 16 values of average IRS's are 
called Average Return Samples (ARS) and are computed using an exponentially decaying 
averaging technique over a period of approximately 2 seconds. It is the slope of the 
ARS's which is examined to determine significant wave height. 

When a square pulse impinges upon an ocean surface with non-negligible significant 
wave height, the shape of the mean return pulse is altered. The geometry of a sea surface 
with non-negligible significant wave heights is illustrated in figure 10. Nate that the 
crests of the waves are illuminated prior to the time at which the calm sea would have 
been illuminated. Similarly, the troughs of the waves are not illuminated until after the 
time at which the calm sea would have been illuminated. The net result of these effects 
is that the mean power received for non-negligible sea state does not reach its foil pla- 
teau value until after the time at which the mean power received from a calm sea reaches 
its plateau values. The leading edge of the mean pulse shape for several non-negligible 
sea states is characterized in figure 11. From this figure it can be seen that the slope 
of the mean return pulse is related to significant wave height. If the mean pulse shape 
for negligible sea state was known precisely, then the significant wave height could be 
determined by analyzing the departure of the mean return pulse for the non-negligible sea 
state from the mean return pulse shape for calm sea. 

It has been shown (refs. 15, 16, 17) that the mean return waveform can be treated 
as a convolution of the composite of the transmitted pulse and transmitter and receiv- 
er bandwidth effects, the non-coherent surface (calm sea) impulse response, the ocean 
surface height probability density function, and the tracking loop jitter. The bandwidth 
effects and the probability density function can be modeled as Gaussian distributions, 
the calm sea impulse can be modeled as a step function, and the tracking loop jitter is 
assumed to be unaffected by sea state. 
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Assuming that pointing angle errors have negligible effect: upon the leading edge of 
the waveform, Brown and Miller (ref. IB) have shown that a gcxci appioximation of the re- 
turn power can be represented as a Gaussian function 


y(t) - a P (~) + d 


( 15 ) 


where a is the return waveform amplitude; b is the time origin; c is the return wave- 
form risetime; d is the return waveform baseline amplitude. The evaluation of the func- 
tion P in equation 15 is performed by: 



~ exp (-q 2 /2) dg 


(16) 


Therefore , it is necessary bo estimate the four parameters a, b, c, and d from which the 
slope of the return waveform, and in turn the significant wave height , can be determined. 
The technique used to estimate the four-parameter function y(t) is the method of least 
separes (ref. 19). 

The convolution model for the mean return waveform assumes that the risetime param- 
eter c is a composite of two Gaussian distributions: the calm sea pulsewidth a ■ and the 

rough sea rms a , both expressed in nanoseconds. Since the convolution of two Gaussian 
distributions is itself a Gaussian distribution, 


2 2 

a + a 
s c 


(17) 


The significant wave height is defined as four times the rms sea height relative to mean 
sea level. Conversion to units of meters by multiplying by the two-way speed of light 
yields 


Sm = 4a g = 0.6 



(18) 


Significant wave heights have been computed on thousands of passes of GEOS- 3 data 
and it has been found that the algorithm typically converges within 2 or 3 iterations. In 
addition, it has been determined that the final converged estimate of significant wave 
height is not particularly sensitive to the initial guesses of a Q , b Q , c , and d Q . In 
practice, the same initial guesses are used for all passes. Eor each significant wave 
height estimate after the first frame of data, the converged values of a, b, c, and d for 
the previous frame of data are used as the initial guesses for a b , c , and d . There- 
fore, it is necessary to have an accurate estimate of the calm sea pulse width (a ) in 
order to calculate significant wave height using equation (18). Early in the GEOS-3 
mission, many passes of the satellite over areas where ship measurements indicated calm 
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seas were analyzed in order to determine that value of which would yield an estimate of 
SWH = 0 for those passes. The value arrival at was 7.49 nanoseconds. This value has been 
examined (refs. 9, 10) a*v- found to be accurate. 

Due to the electronic characteristics of the altimeter, the mean return waveform 
represented by the ARS's contains noise, waveforms from two adjacent data frames can 
differ substantially, even though the altimeter is receiving data from ocean areas only 
a few kilometers apart. This feet, combined with the inherent errors introduced in the 
estimation process, can produce an estimate of c which causes the term under the radical 
in equation (18) to become negative. This is especially true when the sea state is very 
calm. It should be noted that for moderate values of significant wave height, small 
errors in the estimate of c do not cause large errors in the calculated value of signifi- 
cant wave height. However, the estimate of significant wave height is very sensitive to 
even small errors in c for calm seas. 

When the noise characteristics of the altimeter, the estimation errors, and the 
algebraic sensitivity of the estimate to c are combined, it is obvious that smoothing the 
estimate should provide more confidence in its accuracy, especially when the sea state is 
calm. Empirical studies indicate that the longest segment over which the return pulse can 
be assumed to be correlated is approximately 140 kilometers (or 21 seconds). A sliding 
21-second rectangular filter has been employed by the significant wave height estimation 
software to smooth the estimate. Either the estimate of c or the calculated value of 
significant wave height can be smoothed and the results have been shown (ref. 14) to be 
nearly identical. . Fbr computational ease, the estimate of c is smoothed in this investi- 
gation. Even when the estimate is smoothed, the term under the radical in equation (18) 
can still become negative. Such cases have no physical meaning and a value of SWH = 0 is 
assumed. 


Wind Speed Estimation Algorithm 

A widely used parameter in the analysis of radar returns is o q , which is the radar 
backscattering per unit scattering area and is a measure of the surface as a radar 
reflector. It has been shown (ref. 20) that a Q (0 o ), which is the value of cr for no off- 
nadir pointing angle, may be written as 

.„(»•) - . Jupi- ti9) 

~ Clr 

where is a constant of proportionality which depends upon the density of the slopes of 

2 

the filtered surface, R(0°) is the Fresnel power reflection coefficient for the sea at 

2 

normal incidence and 13,9 GHz, and ^ is the mean square slope of the filtered surface, 
which can be logarithmically related to the wind speed W by 


q T = A ln(W) + B 


( 20 ) 


Assuming that a remains essentially unchanged with variations in wind speed, 

. * o ( 0°) 2 

°o (0 > = A ln(W) + B ( 21 ) 

2 

Rearranging terms and using -2.1 dB for R q (0°) gives 

W = exp (- 1 - 0 \ --- - ) (22) 

a (0° ) + 2.1 

where x = — 2 and a Q (0°) is given in dB. The linear coefficients A and B were 

de ter mined by fitting a curve to the GEOS- 3 a Q (0°) estimates where wind speed ground truth 
was available. These data, shown in figure 12, suggest the use of a cusped curve. The 
least squares determination of these constants (ref. 21) yields different values for wind 
speeds above and below 9.2 m/sec. For wind speeds below 9.2 m/sec, A is 0.02098 and B is 
0.10075, whereas for wind speeds above 9.2 m/sec, A is 0.08289 and B is -0.012664. 

Thus, using equation (22), wind speed can be determined solely as a function of 
a o (0°). It should be noted, however, that the received pcwer decreases in a predictable 
manner as the satellite pointing angle increases from nadir. Since the GEOS- 3 pointing 
angle varies by as much as 0.8°, this effect must be accounted for in the calculation of 
a Q . Brown and Curry (ref. 22) developed a scheme for correcting a Q for pointing angle 
variations and that scheme has been implemented in the wind estimation algorithm. 

RESULTS 

The results of this investigation consist of a series of sea surface topographic 
maps and two sets of histograms, one depicting the Significant Wave Height and the other 
the ocean surface wind speed. 


Altimeter Dat i Set 

For the purpose of this investigation, it has been necessary to select GEOS- 3 
altimeter passes whose ground tracks enter the ELM coverage area, which comprises, for the 
SWH and wind speed study, an area from latitude 27° to 35° N and longitude frcm the U.S. 
coastline to 284° E; and 26° to 37° latitude and a longitude of 2° farther east for sea 
surface topograpliies. 

After the altimeter passes which enter the BLM areas are identified, it is necessary 
to ascertain the altimeter status and the telemetry mode of the data. Any passes when the 
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altimeter status indicates th?t the altimeter is not locked in the tracking mode, usually 
during and shortly after the time when the satellite passes over land, cannot be included 
in the data set. The telemetry mode can be: global, intensive with no waveforms, inten- 

sive with 8 waveforms or intensive with 16 waveforms. If the altimeter is in the global 
mode an accurat sea surface height cannot be estimated for this area due to the noise of 
the measurements. Also no waveform information is reported so that the SWH cannot be 
estimated. If the altimeter is operating in the intensive mode, waveform information may 
be or may not be reported. Since all 16 AES's are required in order to estimate SWH, only 
the intensive mode with 16 waveforms can be accepted for the computation of the SWH. All 
intensive modes can be used to obtain the sea surface heights and all modes can be used to 
calculate the wind speed. The data products described below have been reduced when the 
appropriate data rate and telemetry modes were in effect. 

Sea Surface Topographic Maps 

'r\e end product of the sea surface height measurements obtained by the GEOS- 3 radar 
altimeter and processed as described above is a series of sea surface topographic maps 
which iric ude not only height contours but also graduated arrows indicating the magnitude 
and direction of the surface geos trophic flow* At the same time that the monthly sea 
surface maps are produced, a map shewing the sub-satellite tracks for the passes included 
in each month is prepared and they are shown as figures 13, 14, 15. A three-year mean 
sea surface topographic nap was also produced by averaging six half-year mean reticulated 
data tapes. The average was done in a way similar to the monthly maps. All measurements 
falling into a 1/2 degree cell were averaged and empty cells were interpolated. An indi- 
cation of the total amount of passes in the three years of data collection is found in 
figure 16. Maps of mean sea surface topography created for each month in the three years 
of this study from July 1975 to June 1978 in addition to the above mentioned three year 
mean map appear as Appendix A. 

Significant Wave Height and Wind Speed Histograms 

The histograms for significant wave height and wind speed are presented in Appendices 
B and C. In each figure, the axis of abscissa is graduated in meters and the axis of 
ordinate is graduated in percent. In addition, each bin of the histogram is labeled with 
the actual number of samples included in the bin and each figure contains the number of 
points in the total distribution as well as the mean arri standard deviation of the dis- 
tribution. 

The HLM study region was divided into five areas (see figure 17). Each area 
represents a different sector of the ocean, the areas contain the following environ- 
ments: (1) coastal waters, northern sector; (2) coastal waters, southern sector; (3) 
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Gulf Stream waters, northern sector; (4) Qilf Stream waters, southern sector; and (5) 
Sargasso waters. Notice that the division between northern and southern sectors occurs at 
about the area of the Blake Escarpment where the Gulf Stream changes course, from almost 
due north to more northeasterly. The histograms for each of the five BLM areas and a 
composite histogram of all the areas are given for each month, composite of all years; for 
each season, also a composite of all years; and for the entire study time. For seasonal 
histograms, winter is defined as December, January, and February; spring as March, April, 
and May; summer as June, July, anti August; and fall as September, October, and November. 

DISCUSSION 

The results obtained in this investigation are discussed here. Fbr a more conven- 
ient interpretation, they have been divided into two sub- sections: Sea Surface Topograph- 

ic Maps and Significant Wjvc Height and Wind Speed Histograms. 

Sea Surface Topographic Maps 

TWo conflicting problems immediately appear when attempting to understand the mean 
topographic sea surfaces obtained frcm the individual profiles in the GKOS- 3 data base. 

The first is the number and distribution of passes in the geographic area and the second 
is the time scale resolution of the oceanographic phenomena of interest. Since the GEOS- 3 
radar altimeter did not acquire data continuously over the region of interest because of 
scheduling and power constraints, the number and distribution of passes in the Qilf Stream 
area a ic somewhat random during mast time periods. Ibis constraint, forces the number of 
passes available per unit of time to be highly variable. Figure 18 shows that the distri- 
bution of crossing points ranges frcm 0 to 76 per month. With no crossing points, any 
analysis using the data from January, 1976, will be futile. The data from June 1977, on 
the other hand, shows much premise because with 28 crossings, a reasonable monthly topo- 
graphic surface is computed. Using the data depicted for April 1977, with the densest 
coverage for any month, 76 crossing points, results in a reliable surface. Therefore, 
before estimating any surface over a specified time interval, the distribution statistics 
and the distribution and number of crossing points of the passes should be evaluated prior 
to making interpretations of the results. 

The following is a series of short explanatory comments on the sea surface topograph- 
ic maps produced for this study: 

July 1975 There is a well-distributed coverage of north-to-south passes integrated 
with a dense coverage of south- to- north passes in the northern part, of the study area and 
one pass in the southern reaches. A conspi''uous Gulf Stream can be traced throughout the 
entire bight with a well defined Charleston Bump. A marked depression at 33° N, 284° E 
could be an tody soon to be absorbed by the main current. 
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Auvjust 107*' Hiore is a fairly dense average n tto ivntral puts of tto bight with 
an untied piss in the southern reaches. Hie ext rom • nortlH'rn pirts arc thinly a vend by 
a ixirth-to-south piss. Hie Gulf Stream can bo tract*! in tin* entire bight, However, litt le 
emphasis slxxild bo placed on t ho southern pirt where tho unt iod piss may have contributed 
to a veiy slurp slop'. Hie ami around tto Charleston buni> is wt'll do final by a small 
hump. 

Kept qjixvr 1078 The coverage for this month is wt'll distrihdai with almost no ami 
lx'ttor tlvm another. Hie main Ixxly of tto Calf Stream can to traaxl thnxighout the bight , 
losing jdnrpiess tuwanl the extreme northern roadies. Hie Charleston Pump atipears as a 
flattish lunate area. Hie most distinguisliable feature is a 80 cm depression located at 
30° N, 287° R wti ! ch lias toon identified as a cold water eddy. 

(Vtobor 1078 average for this month is adequate; there is a large gap between lxirtli- 
to-south pisses yet this gap is well covernl by soutli-ttv-north passes. Hie only open area 
where caut ion slxxild bo applied appears to to at about 30° N, 281° R wtiere there is a 
largo ami with no avoixigo. Hie Gulf Stream is wall do final except in tlx:' nor t tom 
xxmcIx's where it appears to flow inlaixl. Hiis is probably due to the appearance of a 
baly of cold wuter storewatxl of the Gulf Stream. A mid water aldy can to idont i fiat at 
20.8° N, 282° R. 

Nove mtor 1073 Bassos during this month form a very well distributed grid with no 
distinguishable portions of poor coverage. Hie Gulf Stream is well definal with a sur- 
prising absence of the Charleston Bump. Hie cold Lxxly at Onslow Bay is still present. 

toain tor 1078 Nortli- to- south coverage is goal as is the scxitli-to-north although they 
do not combine to pmlua' a symmetrical grid. Hiis month is marked by fxxir zones of veiy 
add voter storowanl of the Gulf Stream aixl a pronounced storowanl shift of tto Gulf 
St -reum mxir Savannah. 

January 107 h Hiere is a fairly even distribution of nortlv-to-south passes; unfortu- 
nately*, there wore no soutli-to-north pisses to tie it down. Analysis from this month 
slxxild to appixudxxi cunt iously at Ix'St. Due to those problems, there is an evident 
northoast-soutliwost elongation and a protiouiv.xxl depression at alxiut - 30° N, 282.8° R. Hiis 
pis it. ion is auncident with the path of an aldy, but it could also Lx' mincidental, 

Februa ry lQ 7 _ (i Altlxxigh only six passes are obtained, the spicing of tto north- to- 
saith pisses aixl tiie fact that these are tied by two soutli-to-north passes Ix'lp to obtain 
a somewhat clearer topxiraphy. Nonetheless, caution is again exercisal. The Gulf Sr earn 
is fairly evident ltd then' are no eddies present. 

torch 107 (1 This is another px'r month regarding the number aixl li vat ion of pisses. 

Bart iailarly pxir in average aixl immediately evident in the topographic nvip is tiie ox- 
t mne saithern pirti>- where there is a lack of data between pisses. 

April lQ7 t> Here again, altlxxigh tto numtor of pisses is mit large, their spicing 
is goal ontxigh to piatuce a fairly' accurate map. Hie am them sector again stows signs 
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of a systematic interpolation. The Gulf Stream shows fairly well with a pronounce! 
Charleston Bump. Oily one eddy is evident, located at about 30.5° N, 282° E. 

May 1976 This month appears to be well covered by both north-to-south and south-to- 
north passes. Unfortunately, the most shoreward of the north-to-south passes is not tied 
to the south- to-north passes thus leaving a fairly large open space close to shore. No 
eddies are apparent, ard the Gulf Stream is observed throughout, the bight. 

Ju ne 1976 Although not too many passes cover the area, those present, are fairly well 
distributed producing a good coverage. The topographic map shows many mare features than 
the previous several months. Hie Gulf Stream is traceable up to about 33° N where it 
becomes diffused. A large oold water depression is evident at 33° N, 282.5° E and an 
eddy is observed at 30.5° N, 281.5° E. 

July 1976 There is a good distribution of passes with only one south-to-north pass 
untied to the rest, of the system. The topographic map shows a clearly defined Gulf Stream, 
twr cold water areas close to shore, and two eddies. The Charleston Bump is well defined 
and a smaller bump is also evidenced at 29.5° N, 281° E. 

August 1976 A very good spacing of passes is evidenced this month. TWo south-to- 
north passes in the extreme southern reaches are not tied to the rest. The topographic 
map shows a well defined Gulf Stream running the entire length of the map with a strong 
showing of the Charleston Bump. There are no eddies visible. 

September 1976 The passes within this month are well distributed, and only in the 
extreme southern reaches is there an untied pass. The Gulf Stream is well defined 
throughout the map and the Charleston Bump is well marked. 

October 1976 Coverage for this month is more than adequate. The only less desirable 
feature is a large void in the north-to-south passes fairly close to land. Yet this 
deficiency is more than made up by a good number of crossing south-to-north passes. The 
topographic map shows a well marked Gulf Stream and a very cold depression shoreward at 
about long Bay. 

November 1976 Coverage for this month is only fair. Spacing between north-to-south 
passes is too wide except close to shore. South-to-north passes are absent in the north- 
ern reaches. The Gulf Stream is visible throughout and the depression off Long Bay is 
still shown. 

December 1976 The passes in this month are almost identical to those in last month 
providing only a fair coverage. The Gulf Stream is observed fairly well. There is a oon- 
picuous absence of the Charleston Bump and the depression off Long Bay is still present. 

January 1977 Although the number of passes is not large, the spacing is fairly 
symmetrical providing a very good coverage. The Gulf Stream is very distinct throughout. 
The depression off Long Bay is still shewn and no eddies are identified. 

F^truary 1977 Very good coverage is again provided for this month and only one 
south-to-north pass near the extreme end is untied. The Gulf Stream is again shown 
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throughout. The depression at Long Bay is shifted north and the Charleston Bump is very 
conspicuous. 

March 1977 Again, coverage for this month is cfuite good with two untied passes in 
the southern reaches. Hie Gulf Stream is readily observed. Hie depression at long Bay 
seen in the last few months has new moved northward into Onslow Bay. A noticeable 
Charleston Rump is shown immediately west of a cold water eddy located at 30.5° N, 

281.5° E. 

April 1977 Coverage for this month is excellent, due not only to a large number of 
passes but also to the spacing of these passes. Hie Gulf Stream is well defined with a 
very interesting loop at about 30° N, 279.5° E. Hie depression off Onslow Bay has moved 
fartlier north and has intensified. Hie Charleston Bump is shewn much farther out at sea 
than ever observed. It is shown at about 30° N, 282° E rather than at a more customary 
30.5° N, 281° E. 

May 1977 Large spaces between south- to-north passes reduce the giality of this 
month* s data. Hie Gulf Stream is not readily followed in its entirety due to a large 
Charleston Bump which could be due to the lack of passes in that area. Hie depression 
off Cnslow Bay has again moved south to the proximity of Long Bay. 

June 1977 Hie month of June is also of poor quality in coverage due to the poor 
spacing of the passes. The Charleston Bump appears as a large elongated feature extend- 
ing well offshore. The depression off Long Bay has now intensified. 

July 1977 Good coverage for this month is provided by an adequate supply of south- 
to-north passes crossing the north- to- south passes which lack a pass through the center 
of the area. Hie topographic map shows a well developed Gulf Stream, the Charleston 
Bump and the Long Bay depression. No eddies are visible. 

August 1977 A good spacing of passes provides very good coverage for this month. 

Hie Gulf Stream is well described and so is the Charleston Bump. A noticeable feature is 
what appears to be a cold oddy just east of the Bump and located at about 30.5° N, 281° E. 
This depression is about 60 cm deep. 

September 1977 Coverage for this month is only fair, mainly because there are only 
three south- to-north passes and these are somewhat north of center. The topographic nap 
shews a well developed Gulf Stream but it is not well described south of about 30° since 
there is no coverage there; thus, the values are interpolated. 

October 1977 Passes within this month are well spaced; unfortunately, some of the 
south- to-north passes stop short of the coastline. Hie Gulf Stream for this month pre- 
sents an interesting anomaly at 30° N, 280.5° E. This anomaly, a bump of more than 40 cm 
in height, is well off the location of the Charleston Bump. 

November 1977 Coverage for this month is adequate, but there is a large spacing 
between south- to-north passes. The topographic map reflects this void; the Gulf Stream 
is not continuous. 
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December 1977 The passes comprising this month are not only few, but their spacing 
yields only nine crossing points. The topographic map reflects this quality. In the 
southern reaches where there are six of the nine crossing points, the Gulf Stream is well 
defined. In the northern reaches, it is diffused. 

January 1978 The lack of passes for this month, and the arrangement of these passes 
which yield only three crossings, result in a very poor topographic map. Oily the very 
general circulation pattern can be identified. 

February 1978 This is again, a poor month for coverage. Attention should be placed 
at the close contour at the extreme southern parts. These are the product of automatic 
interpolation due to a lack of data. 

March 1978 The number of passes available for this month produces a slightly better 
coverage. Eue to the lack of passes in the southern sections, the Gulf Stream is not 
modelled south of 32°. The most noticeable feature is the large depression off Long Bay. 

April 1978 A good coverage is provided by a fairly dense spacing of north- to- south 
passes with only four north- to- south passes -tying the system. The topographic map shows 
a Gulf Stream which is traceable throughout. A large depression off Long Bay is now 
accompanied by another depression off Onslow Ba,. 

May 1978 This is again a poorly covered month with poor quality due to the arrange- 
ment of passes. There is a very wide open space in the center of the region. The result 
of this coverage is readily apparent in the topographic map. The almost north-south 
contour lines in the center of the map should be interpreted with extreme caution as they 
may be the product of automatically interpolating values where there are no measurements. 

June 1978 This month presents a good coverage as passes in both directions are well 
distributed. The topographic map shews a distinct Gulf Stream. A Charleston Bump is 
displaced south and the cold water depressions off Long and Onslow Bays are present. 

The resulting three-year mean surface contains several interesting oceanographic 
and geodetic areas. First, the Gulf Stream is very well-defined from the Florida Straits 
to approximately 33° N including the effects due to the Charleston Bump just below 32° N. 
The lack of a streamline just east of Cape Hatteras is rather puzzling. TWo explanations 
are possible. First, it is believed that the reference geoid is not able to resolve the 
steep gravity anomaly in that region where the Gulf Stream also tends to hug the shelf 
break. Second, just southeast of Cape Hatteras there seems to be a semi-permanent eddy 
field produced by the coalescing eddies (ref. 23). In a survey of the maps published by 
NOAA in GULFSTREAM covering the same period, eddies were found to exist in this area the 
majority of the time. This may be responsible for the depression southeast of Cape 
Hatteras and it may also cause the distortion of the contours that delineate the Gulf 
Stream. The relatively flat surface in the Sargasso Sea is as expected except for the 
deme centered at 27° N, 284° E. This large done results from three different problems: 
the referenced geoid, the data processing, and a ccnbination of pass distribution and 
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interpolation. First, the geoid has been found to read low by approximately 80 cm compared 
to an altimetric geoid derived by Marsh et al. (ref. 24). Second, sections of the north- 
bound passes in the same area were edited when orienting the individual profiles. However, 
the edited sections were included in the creation of the surface nap thus leading to the 
done effect. And third, the east-west orientation of the southern streamline results from 
missing data in the Bahama Islands region. Those points were then interpolated using pro- 
files on each side of the Islands to complete the gridded surface. 

Significant Whve Height and Wind Speed Histograms 

The monthly significant wave height and wind speed histograms produced during this 
study have been cursorily examined and found to be deficient in a statistical sense for 
characterizing the oceanic conditions on a month-to-month basis during the GEOS- 3 data 
collection period. Ftor the 36 months between July 1975, and June 1978, the number of 
measurements of wind speed per month in the South Atlantic Bight varied from 86 to 801. 

The average number of samples was 360 but there ware 10 months When the number was less 
than 170. In each individual BLM sub-area the number would be of course much smaller. 
Additionally, the telemetry mode settings needed for the collection of SWH measurements 
were not programmed as often as were the settings needed for windspeed measurements . 

Hence, the number of samples available for the monthly histograms would be less for the 
SWH data set. For wind speed analysis, all BLM areas are combined into one which yields, 
for the worst month a total of 86 samples. The credibility of that month's histogram can 
be characterized as follows. Given the inclination of the satellite's orbit, the maximum 
length of a GEOS-3 ground track through the BLM composite area is about 800 km. For a 
sampling period of 3.2 seconds and an effective groundspeed of approximately 7 km/sec, it 
is possible that about 86 samples be taken during a single orbit. With only 86 samples 
in a histogram, it is possible that all data were collected from only a few orbits. If 
these occurred during atypical weather conditions then that month's data is not character- 
istic of the entire month's sea conditions. 

Several of the monthly data sets are most likely undersampled. Others are undoubtedly 
truly representative of the environmental conditions. To remain conservative, it is appar- 
ent that averaging for more than one month is desirable. But, for maximum impact it is 
desirable that the GEOS-3 data be used to describe the sea state and wind speed in a 
continuous fashion for the entire GEOS-3 lifetime. To accomplish this, and to superfi- 
cially characterize the effects of geographic location within the South Atlantic Bight the 
following three kinds of analysis have been performed in which spatial and temporal aver- 
aging ware employed. 

The average monthly histograms have been studied and they reveal the annual variations 
in significant wave height and wind speed for the composite BLM area. Figure 19 shows 
the variation of the mean and variance for these histograms. Both means are snallest in 
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the simmer and largest in the winter as expected. The distribution mean for wind speed 
varies through a much larger range than does the SWH mean. This will be discussed in more 
detail below. The SWH and wind speed distribution variances seem to undergo similar cycles 
through the year hut a physical explanation of their behavior will be discussed also in the 
description of the second analysis. 

With only about three years of data available lor analysis, it is reasonable to 
assume that the year-to-year variability of the environmental conditions is large enough 
to render the preceding average annual cycle less than optimum. Therefore, to keep the 
sample size of the SWH and wind speed histograms sufficiently large and yet still show the 
year-to-year fluctuations, a sliding three-month window was used in a second kind of 
analysis. The behavior of the means of the resulting SWH and wind speed distributions is 
shown in figure 20. The annual cycles are clearly seen. Although there is a double maxi- 
mum in the SWH variability in the fall and winter of 1978, these cycles are generally 
single-peaked. The small irregularities seen from year to year are, of course, masked in 
the annual average cycle shown in figure 19 . Figure 21 shews the corresponding variances 
for the sliding window distributions. The SWH histogram variance does not undergo well- 
defined yearly cycles, perhaps because the range of sea states common to the South Atlantic 
Bight is small. The distribution for mid-winter may be only 30-40 cm wider than is the 
distribution for mid-summer. This, in turn, is only a fraction of the width of the bins 
used in these analyses. For the wind speed data, however, the situation is much clearer. 
The variance is minimal in mid-summer, but there appears to be two maxima for each year's 
data record. These are centered around the maxima of the distribution means shown in 
figure 20. This can be interpreted in the following fashion: when the large-scale weath- 

er patterns change in the fall, there are periods of time in which the environmental con- 
ditions still resemble those of the sumxer. In between, hewever, storm systems begin to 
move through the area increasing the wind speed which, if substained for a sufficient 
period of time, also increases the SWH. Hence, the range of SWH and wind speed values 
occurring in this season is increased reflecting these two different weather regimes. In 
terms of the wind speed histogram, at least, this appears as an increase in the variance. 

By mid-winter, there are fewer and fewer mild periods so that although the mean wind speed 
continues to increase to its naximum, the variance drops. In springtime, the process is 
repeated in the reverse direction. 

The final kind of analysis pertains to the spatial variability of SWH and wind speed. 
To avoid undersampling, it is decided that each BLM area can be studied if sufficient 
temporal averaging is done. Hence, the monthly histograms for each area are added to 
those from other years, and the means of the resulting average monthly histograms are 
plotted in figure 22 for the wind speed data. For SWH, no detectable difference in dis- 
tribution means or variances are found between the BLM areas. For wind speed there is an 
apparent difference. The mean wind speed for area 3 consistently exceeded the means for 
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the other areas. No explanation will be tendered but the well-known prevalence of high 
winds near Cape Hatteras, North Carolina suggests that the meteorological conditions of 
the region are accurately reflected in these histogram analyses. Hie GEOS- 3 data set used 
in this study can easily be analyzed in other ways and in various degrees of detail. The 
analyses above are merely examples of the kind of climatological studies that are possible. 

CONCLUSIONS 

Hie analysis of long-term altimeter data for sea surface topographies leads to the 
following conclusions. Due to the various uncertainties in geoid, orbit, and intrinsic 
noises, the mean surface topographic naps are to be evaluated with extreme caution. Among 
the uncertainties, the accuracy of the geoid is *.ie most critical element because the geoid 
anomalies can have, not only the same wavelength as ocean dynamic phenomena, but the geoid 
signal magnitude can also be rruch higher. Hie difference between the long-term (three 
years) and any short-term mean (one month) reveals the changes of Gulf Stream path and 
some eddy movements. In principle, the uncertainty of the geoid should not be a problem 
here, but it still causes difficulty. In this case, difficulty comes from the data density 
from the altimeter. Hie long-term mean is produced by using the average of all the ad- 
justed heights within a 1/2 degree square. Hiere could be considerable geoid change with- 
in such an area. Hie averaging process reduces the true variations. 

Furthermore, it is apparent that seasonal and annual variations in significant wave 
height and wind speed are readily detectable using waveform data from satellite radar 
altimetry. Hiis is demonstrated with GEOS- 3 measurements over the South Atlantic Bight 
for this study. Different degrees of spatial and temporal averaging are used to ensure 
that under sampling is avoided. On the smallest spatial scale tested, seme differences 
in the resultant wind speed distributions are found to agree with the known climatology 
of the area. No wave height spatial difference is detected, an expected observation 
given the low wave heights found in the South Atlantic Bight, 
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Figure 3. Overall Sea Surface Height Data Processing 
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GULF STREAM ANALYSIS DATA PROCESSING 



Preliminary Data Processing by DOS Computer Program 
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Figure 6. Sub-band Analysis for the Determination of the Position of the 
Termination of the Straight Line Fit 




Figure 7. Comparison of Bill and Thinnei-out Sea Surface Profile 
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Figure 15* GEOS- 3 Sub-satellite Tracks of Passes Used for the Monthly Sea 
Surface Topography July 1977 to June 1978 
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Figure 17. RLM SWH and Wind Speed Investigation Areas 
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Figure 18. Monthly Frequency Distribution of the Nurifcer 
GEOS- 3 Crossings Over the Study Area 
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Figure 19. Means and Variances for the Average Monthly Significant ft&ve Height (SWH) 
and Mean Wind Speed (W) Distributions for the GEOS- 3 Data Set 
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Figure 21. 
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MEAN - 8. 2 > SIGMA « 2. 8 
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DISTRIBUTION OF WIND SPEED IN METERS PER SECOND 

SPRING SEASON 
BLM AREA #1 

N - 344. MEAN ■ 5. 2, SIGMA * 3. 5 



DISTRIBUTION OF WIND SPEED IN METERS PER SECOND 

SPRING SEASON 
BLM AREA #2 

N = 299. MEAN = 4. 5, SIGMA = 2. & 
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DISTRIBUTION OF 

WIND SPEED 

IN METERS PER SECOND 


SPRING SEASON 


PLM AREA 

#3 

N » 438. 

MEAN ■ 7. 3, 

SIGMA >3.4 


DISTRIBUTION OF WIND SPEED IN METERS PER SECOND 

SPRING SEASON 
BLM AREA «4 

N - 492, MEAN - 6. 2, SIGMA - 3. 2 
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DISTRIBUTION OF WIND SPEED IN METERS PER SECOND 

SPRING SEASON 
BLM AREA #5 

N - 1237, MEAN - 6. 3. SIGMA - 3. 2 
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DISTRIBUTION OF WIND SPEED IN METERS PER SECOND 

SPRING SEASON 
ALL BLM AREAS 

N = 2770, MEAN ■ 6. 1. SIGMA * 3. 3 
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DISTRIBUTION OF WIND SPEED IN METERS PER SECOND 

SUMMER SEASON 
BLM AREA #5 

N = 2499, MEAN « 5. 2. SIGMA = 3. 0 



DISTRIBUTION OF WIND SPEED IN METERS PER SECOND 

SUMMER SEASON 
ALL BLM AREAS 

N = 5492, MEAN * 5. 1, SIGMA = 2. 9 
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DISTRIBUTION OF MIND SPEED IN METERS PER SECOND 

FALL SEASON 
BLM AREA *1 

N « 627. MEAN = 5. 9, SIGMA - 3. 5 
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DISTRIBUTION OF WIND SPEED IN METERS PER SECOND 

FALL SEASON 
BLM AREA *3 

N « 783, MEAN * 6. 7, SIGMA = 3. 5 



DISTRIBUTION OF WIND SPEED IN METERS PER SECOND 

FALL SEASON 
BLM AREA *4 

N = 810, MEAN = 5. 9, SIGMA = 3. 2 
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DISTRIBUTION OF WIND SPEED IN METERS PER SECOND 

FALL SEASON 
BLM AREA #5 

N - 2351, MEAN * 5. 8. SIGMA = 3. 3 
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DISTRIBUTION OF WIND SPEED IN METERS PER SECOND 

FALL SEASON 
ALL BLM AREAS 

N = 5060 . MEAN = 6. 0, SIGMA = 3. 4 
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DISTRIBUTION OF WIND SPEED IN METERS PER SECOND 

WINTER SEASON 
BLM AREA #3 

N *> 352, MEAN « 9. 0, SIGMA ® 2. 6 



DISTRIBUTION OF WIND SPEED IN METERS PER SECOND 
WINTER SEASON 
BLM AREA «4 

N = 374, MEAN ® 8. 8, SIGMA - 2. 5 
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DISTRIBUTION OF WIND SPEED IN METERS PER SECOND 
ENTIRE MISSION 
BLM AREA #4 

N = 2537, MEAN = 6. 0, SIGMA = 3. 3 
C 50 
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DISTRIBUTION OF WIND SPEED IN METERS PER SECOND 
ENTIRE MISSION 
BLM AREA «5 

N = 7090, MEAN = 6. 0, SIGMA = 3. 4 
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